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Abstract 


Raindrops measured by two orthogonal cameras were classified by shape and orientation 
to determine the nature of the oscillation. A physical model based on potential energy was 
then developed to study the amplitude variation of oscillating drops. The model results 
show that oscillations occur about the equilibrium axis ratio, but the time average axis 
ratio is significantly more spherical for large amplitudes because of asymmetry In the 
surface potential energy. A generalization of the model to oscillations produced by turbu- 
lence yields average axis ratios that are consistent with the camera measurements. The 
mode] results for average axis ratios have been applied to rainfall studies with a dual- 
polarized ( vertlcal/horlzontal) radar. 


Introduction 


The equilibrium shape of a falling drop Is the result of a balance between surface 
tension, aerodynamic forces, and hydrostatic and internal pressures. Drops which have been 
observed In laminar flow wind tunnels ‘ and in still air'^ assume shapes that are consistent 
with the equilibrium calculations of Pruppacher and Fitter^. Under turbulent conditions, 
however, drops are not quiescent since they undergo continuous oscillations'* > ^ . In fact, 
photographs of raindrops In the atmosphere Indicate that the average shape departs signifi- 
cantly from equilibrium®. Although these departures ha\e usually been neglected, they take 
on a renewed Importance with recent advances In radar mjteorology that exploit the shape 
and orientation of scatterers to obtain a quantitative or qualitative description of 
precipitation^ . 

The following study consisted of two major efforts. The first was to determine raindrop 
shape from the photographic record, and the second was to construct a physical model of the 
changing shape to compute the time average axis ratio for the analysis of dual-polarization 
radar data. The application of the model is discussed in another paper®. 

Data analysis 

In the late 1950's Jones studied the shape of natural falling raindrops using two cameras 
with optical axes 90° apart®. From photographs he was able to determine the axial ratio and 
volume of 1783 drop.?. Jones' data has been replotted in Figure 1 to show the average axiai 
ratio (vertlcal/horlzontal) with the 955! confidence interval as a function of equivalent 
volume diameter. The solid line In Figure 1 represents the equilibrium drop shape cal- 
culated by Pruppacher and Fitter®. As the 955! confidence Intervals shows, the mean axial 
ratio for raindrops is significantly more spherical than equilibrium. Raindrops larger 
than 4 mm. were also found, to be more spherical than equilibrium, but were not photographed 
In adequate numbers to yield a good estimate of the mean. Although the total sample of I 783 
raindrops Is relatively modest, a preliminary inspection of measurements using a single 
camera’ with several hundred thousand raindrops between 2 and 6 mm diameter shows similar 
results . 

Data from Jones was re-analyzed to determine the nature of raindrop distortions. A 
simple ellipsoidal shape was assumed which is consistent with an oscillation of the funda- 
mental mode. Rotation was not considered as a cause of drop distortion since the energy 
required is much larger than for oscillations'®. Each data point was analyzed for charac- 
teristic shape by comparing the measured vertical and horizontal dimensions. Smaller drops 
were found to be predominately spheres and larger drops predominant i v ellipsoids. It was 
apparent from the shape analysis that raindrops do not oscillate Just in the vertical axlsym 
metric mode because of the significant percentage of horizontal ellipsoids. 

The type of oscillation was examined by combining shapes suggestive of the vertical 
axlsymmetrlc oscillations and those shapes suggestive of horizontal oscillations. The 
result was consistent with the wind tunnel observations of Nelson and Gokale'' that small 
drops (D $ 3 mm) prefer the axlsymmetrlc oscillation (m » 0 degeneracy of the fundamental 
mode) and of Blanchard’ that large drops (D ^ 5 mm) prefer the horizontal oscillation (m»2) 
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The latter Is apparently enhanced by aerodynamics since the external pressure should 
provide a positive feedback to the horizontal deformation. Support for an aerodynamic 
effect Is found In the observation of stable horizontal "prolates" for very large drops*. 

The drop was assumed to have only surface, gravitational and kinetic energies. The 
minimum In the surface plus gravitational energy was calculated numerically. Agreement 
with wind tunnel data and the numerical result of Pruppacher and Fitter for the equilibrium 
axis ratio was within 1!5 for diameters less than 5 mm (Figure 1) demonsv.ratlng that the 
aerodynamic and hydrodynamic pressures are much less Important than the hydrostatic pres- 
sure In determining the axis ratio. 





The surface plus gravitational energy 
function was used as a potential energy "well" 
for anharmonlc oscillations of an ellipsoidal 
drop. The equation of axlsymmetrlc motion for 
the fundamental m'de was obtained from the 


energy relation 

PEq - PE = KE = m(Va2 + + V<j2)/5 (1) 

by solving foi- the temporal change In axis 
ratio 

AR = 6AR^/3 v^/D (2) 

where 

Va = (5KE/m)l/2 ( + 2 Ar2)"1/2 (3) 

and 

AR = c/a = c/b (i|) 

with a constant volume constraint 

D = 8abc. (5) 

An analytical solution for small amplitude 
yielded the Rayleigh oscillation time (tr) 


using the factor of 1/5 In the kinetic energy 
relation. Linear dissipation was added to the 
calculation by assuming KE = B KE where B was 
determined from the result of Lamb’^ for a 
complete cycle. 


Figure 1. Drop axial ratio as a function 
of equivalent spherical volume diameter. 
Circles are for- raindrops (Jones), and 
the curve for equilibrium shaped water 
drops (Pruppacher and Fitter). Con- 
fidence Intervals for the mean have 
been auded to the original raindrop 
data at the 95 % significance level 
based upon the estimated standard de- 
viation and sample size. Triangles 
are present modei for equilibrium 
shape . 

with gravity (g = 980 cm sec"^^ s 0.998 g 
cm2 sec~l and u = 0, yielded values of on 
turning points. 


cm"3 

for 


Numerical Integrations were made for D » 
0.12 cm to compare with the Navler-Stokes 
simulation of Foote^^ for the axlsymmetrlc 
oscillation without gravity with Py, « 1 g 
and a «= 73 dynes cm~fT^ 'me present model 
the large amplitude oscillation of Foote 
(where ARq = 1/1.7) gives the same oscillation 
times (l.e., 1.04 tr for u “ 0.06 cm2 sec~T 
and 1.06 tr for u = 0), and shows similar 
temporal variations In the surface and kinetic 
energies, and in axis ratio with 57 % of the 
Ime spent as a prolate spheroid. The axis 
ratio averaged over one cycle for the above 
two cases, and also for D «= 0 . 3 and 0.6 cm 
cm~2 and o = '^3 dynes cm~l) for both u = 0.01 
ly slightly less than the mean of the undamped 


Asymmetric oscillations of the fundamental mode (ra = 2 degeneracy) Is currentlj under 
investigation. A set of equations similar to those given above is used to solve the 
temporal behavior of the axis ratio. 


Analytical and numerical solutions are more difficult to obtain than the axlsymmetrlc 
case since oscillations occur In the shape of a general ellipsoid. Preliminary analysis 
indicates that the time average axis ratio for the asymmetric oscillation of a raindrop 
can be approximated by the mean of the turning points defined by the potential energy well. 
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Three candidates for the cause of raindrop oscillations were examined: drop Inter- 
actions, turbulence In the air and pumping by wake shedding. A drop Interaction model was 
used to obtain the oscillation energy from a balance between the average collislonal energy 
and viscous dl slpatlon. The time average axis ratio was found to be considerably larger 
than equilibrium for drops with D J 3 mm In heavy rainfall. In contrast, preliminary 
results for wake forcing Indicate that appreciable oscillations occurs only for D * 2 mm. 

For turbulent forcing the oscillation energy was related to turbulence In the Inertial 
subrange using the distance traveled by a raindrop during an oscillation. This gave a 
unique scaling for the oscillation energy with raindrop size and dissipation rate. Mean 
axis I’atlos for the average oscillation energy were calculated from a balance between Input 
of turbulent energy and dissipation by viscosity. The results provided a reasonable fit to 
the mean axis ratios of the camera data and . sultsble envelope for the observed amplitudes. 
It was concluded that atmospheric turbulence and also drop Interactions can produce 
significant Increases In average axis ratio of raindrops resulting In an appreciable change 
In rada backscattei cross sections. 
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